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Magnetic order in orbital models of the iron pnictides 
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We examine the appearance of the experimentally-observed stripe spin- density- wave magnetic 
order in five different orbital models of the iron pnictide parent compounds. A restricted mean- 
field ansatz is used to determine the magnetic phase diagram of each model. Using the random 
phase approximation, we then check this phase diagram by evaluating the static spin susceptibility 
in the paramagnetic state close to the mean- field phase boundaries. The momenta for which the 
susceptibility is peaked indicate in an unbiased way the actual ordering vector of the nearby mean- 
field state. The dominant orbit ally resolved contributions to the spin susceptibility are also examined 
to determine the origin of the magnetic instability. We find that the observed stripe magnetic order 
is possible in four of the models, but it is extremely sensitive to the degree of the nesting between the 
electron and hole Fermi pockets. In the more realistic five-orbital models, this order competes with 
a strong-coupling incommensurate state which appears to be controlled by details of the electronic 
structure below the Fermi energy. We conclude by discussing the implications of our work for the 
origin of the magnetic order in the pnictides. 



PACS numbers: 71.10.Fd, 74.70.Xa, 75.10.Lp 

I. INTRODUCTION 

The origins of high-T c superconductivity is one of the 
most important problems in contemporary condensed 
matter physics. There are now two main material 
classes in which such a state occurs: the cuprates, and 
the recently-discovered iron pnictides P*l The latter have 
therefore been the subject of an intensive research ef- 
fort over the past few years, spurred by the hope to gain 
insight into the mechanisms of the unconventional super- 
conductivity in both families. 

The proximity of antiferromagnetic (AFM) and su- 
perconducting states in the phase diagram of the pnic- 
tides and the cuprates hints that exchange of virtual spin 
fluctuations is responsible for the pairing in both sys- 
tems P™1 Although this already provides strong motiva- 
tion to study the magnetism of the pnictides, the unusual 
magnetic properties are of interest in their own right EEl 
It is now well established that the undoped cuprates 
are AFM Mott insulators P In contrast, there is convinc- 
ing evidence that the iron pnictide parent compounds 
i^FeAsO and AFe2As2 (R and A are rare-earth and 
alkaline-earth elements, respectively) display a metal- 
lic spin-density- wave (SDW) stateP^ the relatively low 
value of the ordered moment!^ and the appearance of 
incommensurate magnetism with dopingp^ the presence 
of several electron and hole Fermi surface above T/v, 
and their significant reconstruction in the AFM state, 
as revealed by angle resolved photoemission spectroscopy 
(ARPES) and quantum oscillation experiments p^^ the 
compounds continue to show metallic transport proper- 
ties below T/v F™^ and spectroscopic measurements sug- 
gest intermediate correlation strengths Q3E2I 

The microscopic description of the AFM state of the 
parent compounds in the cuprates and the pnictides will 
therefore be qualitatively different. While the t-J model 



gives a good account of the former^ there is currently no 
generally-accepted model of the latter. By treating the 
Fe-As planes as a lattice of localized moments, several au- 
thors^"2U have obtained good agreement with neutron- 
scattering data. 25 Although such a picture may be ap- 
propria te fo r the more strongly correlated iron chalco- 
genidesp^it is difficult to reconcile with the experimen- 
tal evidence for a moderately correlated metallic system 
in the 1111 and 122 pnictides. A more promising starting 
point for the pnictides is the prediction of ab initio cal- 
culations that the nesting between the electron and hole 
Fermi pockets derived from the Fe 3d orbitals is responsi- 
ble for the SDW. 26 27 Theoretical work based on this sce- 
nario can be divided into two schools of thought: those 
which argue that only the nesting is imp or t ant J 2 -^"^ and 
those which also attempt to account for the complicated 
orbital structure of the Fermi surfaces pSHSSJ 

The first approach can be dubbed the "excitonic" 
theory, as it is based upon the excitonic instability of 
a semimetal proposed in the 1960sP^ Within this sce- 
nario, the Coulomb attraction between electrons and 
holes about the nested electron and hole Fermi pock- 
ets, respectively, causes the condensation of interband 
electron- hole pairs (excitons). The condensed system is 
a density wave state with ordering vector Q equal to 
the nesting vector of the Fermi surfaces, with secondary 
interactions stabilizing a SDW) 29 * 70 ! Using suitably cho- 
sen phenomenological bands and effective interactions, 
the excitonic theory naturally explains the reconstruc- 
tion of t he Ferm i surface and the continued metalicity 
below Tat j 33 * 37 *^ it can reproduce the results of neutron- 
scattering experiments, 1 34 1 36 1 and it appears to capture 
key aspects of the competition between superconductiv- 
ity and the AFM state l 29 ' 31 ' 33 ' 35 ' 3 ^ 

Despite the success of the excitonic approach, it is un- 
satisfying as a microscopic model of the iron pnictides. 



Specifically, it ignores the complicated orbital structure 
of the different Fermi surfaces j 13 * 71 " ^ and the use of effec- 
tive interactions obscur es the roles of the various intra-ion 
interaction terms! 29 * 33 ' Although the basic physics driving 
the AFM is expected to remain the same, it is desirable to 
include the orbital physics for a number of reasons: the 
variability in superconducting properties of the doped 
compounds suggests that small details of the electronic 
structure could play a significant role in the physicsp' 
the likely presence of moderate correlations implies that 
so me loc al physics should be included in theoretical mod- 
els j 75 * 76 ! while the real-space structure of the orbital wave- 
functions, and the observed orbital reconstruction of the 
Fermi surface below T/v r^ indicates a strong coupling of 
the orbital and magnetic degrees of freedom. 

A large number of different orbital tight- 
binding models have therefore been proposed 
involving between two or five of the Fe 3d or- 
bital S l 4 ^^ l 5i l 52 l 54 l 56 l ^ l a "^ l a5 l ^^ or also including 
the As 4p orbitals 45 80 and even orbitals from outside 
the FeAs planes.^ Although the existence of super- 
conductivity intheseinodels has been extensively 
stud[ed MMMMmMMMMmm the magnetic be- 
haviour remains rather poorly understood. The most 
popular approach to the appearance of AFM order in 
these models is to examine the ground state using a 
standard mean-field an satz that allows for at most two- 
site magnetic unit cells j 47|50 | 52 | 54 | 58 | 59 | 61 | 64 :i very recently, 
this h as be en generalized to a Gutzwiller mean-field 
theoryP"^' As all of the orbital models have a rather 
complicated electronic structure, however, it is by no 
means certain that such phases have the lowest free en- 
ergy, i.e. states with larger unit cells or incommensurate 
ordering vectors cannot be a priori excluded. Indeed, in 
several cases more advanced techniques have been used 
to study the magnetic properties with mixed results: 
whereas in some models the findin gs are c onsistent 
with the usual mean- field ansatz E5 EZE5E1EZ I striking 
divergences have been found in others. 64 The existence 
and robustness of the experimentally-observed stripe 
AFM order in these models at parent compound filling 
therefore remains an important open problem, as it is 
questionable to use them to examine the superconduc- 
tivity of the doped system if they do not give the correct 
magnetic behaviour of the parent compound. 

It is hence desirable to examine the magnetic proper- 
ties of these orbital models in an unbiased way. This can 
be achieved by examining the behaviour of the static spin 
susceptibility Xs(q 5 ^ — 0) m the paramagnetic (PM) 
state: as the temperature is lowered towards the critical 
temperature T mag of a magnetically-ordered phase, the 
susceptibility will diverge at the ordering vector q = Q. 
Practically, we can implement this scheme in a weak- 
coupling approach by first calculating the phase diagram 
of the model using a two-site mean- field ansatz, and 
then examining the location of the highest peaks in the 
random-phase approximation (RPA) spin susceptibility 
at a temperature just above T mag to check the validity of 



the low-temperature mean-field state. Clues to the ori- 
gin of the magnetic instability can be extracted from the 
dominant contributions to the total static susceptibility 
from the orbitally resolved susceptibilities. 

It is impossible to study each of the very many dif- 
ferent orbital models proposed for the pnictides, and so 
in this paper we will restrict ourselves to five models: 
the two-orbital model of Raghu et a/., 43 the three-orbital 
model of Daghofer et a/.p^the four-orbital model of Yu 
et al., bA and the five-orbital models of Kuroki et a/.™ and 
Graser et al..^ The main goal of this paper is to verify 
the appearance and explore the robustness of the stripe 
AFM order in the weak-coupling limit in these very differ- 
ent models, which have been selected as a representative 
sample of the different levels of sophistication available 
in the literature. By this survey we not only hope to 
test the suitability of these models and the weak-coupling 
theory to describe the pnictide parent compounds, but 
also to gain insight into the origin of the AFM state and 
the features of the electronic structure that enhance or 
suppress it. We additionally hope to identify an orbital 
model which gives the required AFM state and is also 
consistent with key results of ab initio calculations and 
ARPES data. Finally, the role of Hartree shifts in renor- 
malizing the band structure and the implications for the 
magnetism is examined. 

Particular attention in our analysis will be paid to the 
five-orbital models of Kuroki et at. and Graser et al. be- 
cause they are currently the most intensively-studied, the 
most sophisticated, and likely most accurate as they both 
result from fits to ab initio band structures of LaFeAsO. 
We note that other five-orbital models of LaFeAsO were 
proposed in Ref. [78] and Ref. [GQl but we will not exam- 
ine them here. The latter is very similar to Kuroki et 
a/.'s model, and therefore we neglect it here in the inter- 
ests of brevity. The former is quite different to Kuroki et 
a/.'s and Graser et a/.'s models as it was not the result 
of fitting to ab initio calculations and does not include 
higher-order hopping terms. This leads to very strong 
ellipticity of the electron Fermi pockets, which substan- 
tially reduces Fermi-surface nesting and should signifi- 
cantly alter the magnetic fluctuations. Since the focus 
of our paper is on nesting and the resulting magnetic 
instabilities, we do not consider this model. 

The structure of the paper is as follows. In section [TT] 
we outline the theoretical basis of our analysis by dis- 
cussing the general form of the Hamiltonian, the mean- 
field decoupling scheme, and the construction of the spin 
susceptibility within RPA. We then proceed to a system- 
atic analysis of the magnetic order in the five models 
In section 



in section III 



IV we compare and contrast 



the ordering properties of the models, and consider the 
implications for the pnictides. We conclude with a sum- 
mary in section W\ 



II. THEORY 



A. Mean-field analysis 



The general Hamiltonian for the orbital models is writ- 
ten H = Hq + Hj. The non-interacting Hamiltonian Hq 
is given by 



fl o = EEE^( k )0 



k,/i,cr J 



(1) 



v^n 



^k^cr 



where dl„„ (d^ uo .) is the creation (annihilation) oper- 



ator for a spin a electron of momentum k in the orbital 
v. T UjfI (k) are the kinetic energy matrix elements for an 
effective tight-binding model defined in t he single- Fe u nit 
cell of the two-dimensional Fe-As planes.l^^^^ 

In all the proposed models, only on-site terms are in- 
cluded in the interaction Hamiltonian 

\ v i u^fi a, a' 

i vj^\i 

+ J 'EE rf t,t4,i rf i,,4 rf i,,,f ( 2 ) 

i i/^/x 

The number and spin operator for the orbital v at 
site i are defined as ri\^^ = d\ va d ivo . and S^ = 

\Ylqq' d[ vc; <;,<;' d ivqf , respectively, where a^^ is the 
vector of the Pauli spin matrices. The matrix elements 
of the kinetic energy of all the models considered here 
were explicitly constructed to obey the orbital and lattice 
symmetries of the Fe-As planes, which has been shown to 
be possible for a one-Fe unit cell in 2-dimensional mod- 
els.^ Even though the interaction terms are purely on- 
site, some restrictions on the parameters still arise from 
the requirement that it should preserve the symmetries of 
the orbital space. The orientation of the axes of the coor- 
dinate system used to define the orbitals must not change 
the symmetries of the Hamiltonian, which implies J = J' 
and V = (2U - 5 J)/4P 

For the sake of brevity we only present results for 
Hund's rule coupling J/U = 0.25. Although this ratio 
is at the upper limit of those employed in the literature, 
it is not unreasonable. Screening can considerably re- 
duce the effective Coulomb interactions U and V from 
their atomic values, which suggests rather large ratios 
of J/U. On the other hand, selecting J < U/3 ensures 
that an electron added to an undoped site pays more en- 
ergy to Coulomb repulsion than it can win from Hund's 
rule, i.e. that the onsite interaction energy suppresses 
charge fluctuations rather than enhancing them. Fur- 
thermore, a mean-field study of a number of models has 
recently found that the ordered moment and the recon- 
structed Fermi surface in the T = OK AFM state is in 
good agreement with experiment for our choice of J/U^ 



When constructing the mean-field phase diagrams of 
the orbital models, we adopt the frequently-employed 

ansatz 54585961 ™ 



( n i,^) 



\e { ^m v 



(3) 



Note that this neglects the possibility of inter-orbital av- 
erages (d\ „ a d- a ) (y 7^ /i), which have been included in 
some studies HZE3 w e restrict ourselves to ferromagnetic 
(FM) and AFM phases with a two-site unit cell, corre- 
sponding to the ordering vectors Q = (0,0), Q = (tt, 0) 
and Q = (7r,7r) in Q, respectively. The Q = (tt,0) AFM 
state gives the magnetic order found in the pnictide par- 
ent compounds. Upon decoupling the interaction term 
Hi, we obtain the mean- field Hamiltonian 

#MF = H o + ^^J2 e v d Kv : * d Kv,v 

k v a 
~ \ E E E ^4+Q,,, A,.,a + NC ,(4) 



k v a 



where N is the number of lattice sites and 



= Un v + (217 - 5 J) J2 ' 



j i± i 



ju#i/ 



Vv 



= Um v + J 2_\ ' 



°fi i 



(5) 
(6) 



H^v 



U 



C = -Uj2nl+-J2ml-(2U-5J)J2n„r 









(7) 



Although the Hartree shifts e v can have rather large 
magnitude, only the relative shifts e v = (5 J — U){n v — 
minjn^}) are important when working at fixed doping. 
These relative shifts are much smaller than the bare 
shifts e v . Note t hat n ot all mean-field studies include 
the Hartree shift sJ=2021 The influence of Hartree shifts on 
the phase diagram will be discussed in section |IVB| 

For given temperature T and interaction constants U 
and J, the stable mean-field state is determined by the 
standard iterative diagonalization technique, where the 
condition of constant charge density 



E< 



1 for 2 orbitals 

2 for 3 and 4 orbitals 

3 for 5 orbitals 



(8) 



is enforced at each iteration step. Upon convergence, the 
free energy per site is calculated. Comparing the free 
energy for each mean-field ansatz allows us to determine 
the stable state of the system. 



B. The magnetic susceptibility 

We define the spin susceptibility by 
X n '(q,ioJ n ) 
= ^E/ dre^ T (r T 5J(q,r)5f(-q,0)), (9) 



V,\l 



where the Fourier-transformed spin operator 5^(q) is re- 
lated to the spin operator at site i by 



The total spin susceptibility is hence written as 

x s (<i,iuj n )= ^2 x JJ (q,^n). (10) 

j=x,y,z 

To obtain the static spin susceptibilities we make the an- 
alytic continuation iuj n — >• uo + z0 + and then take the 
limit cj = 0. Note that it is not necessary to regularize 
the static susceptibility by assuming a non-zero imagi- 
nary part of the frequency Because of the spin-rotation 
invariance of the PM state, we can express Xs(q?^ = 0) 
in terms only of the static transverse susceptibility 

Xs (q^ = 0) = 5 x -+( q ,a; = 0). (11) 

Henceforth we will only be concerned with x h (q> co> = 0) 
as this contains all relevant information on the magnetic 
ordering vector. 

We calculate the transverse spin susceptibility using 
the RPA. Introducing the generalized transverse suscep- 
tibilities 



X„t>u t A<l,iUn) 



1 rf 

k,k' JU 



dre 1 ' 



(12) 
we write the total transverse susceptibility 

x" + (q,^n) = *52xww(<l,iu n ) . (13) 

Summation of the ladder diagrams yields a Dyson equa- 
tion for the generalized susceptibilities 

-+(o) ( ■ \ 

+X^^a,U^^n)K,6,c,dX^,^(q 5 ^n) , (14) 

where the non-zero elements of V a ,b,c,d are 

V a ,a,a,a = U, (15a) 



V a ,a,b,b = J , 




(15b) 


V a ,b,b,a = U - 


-2J, 


(15c) 


V a ,b,a,b = J, 




(15d) 



and we assume a ^ b. The mean-field susceptibilities 
Xvv' u u'faj^n) can be explicitly written as 

-+(o) / • \ 

= - x E E u *y ( k K,/*( k K>' ( k + qK>( k + q) 

k s,s' 

n F (Es y k) -n F (£ s /, k+q ) 



^s,k — ^s',k+q — iVn 



(16) 



where np{E) is the Fermi function, E s ^ are the eigen- 
values of Hmf Q, and the coefficients u s ^ v {k) transform 
the diagonalizing annihilation operators 7 S k correspond- 
ing to the eigenvalues E s ^ into the orbital basis, i.e. 

^,k = E s ^( k h s , k - 

Apart from the total transverse spin susceptibility, 
it is also interesting to consider the dominant contri- 
butions to the sum (13). It is convenient to define 
X~+(q,^ n ) = X^^/Jq^n). These orbitally resolved 
susceptibilities can provide important insight into the re- 
gions of the Fermi surface most strongly involved in the 
magnetic order. We adopt the following enumeration of 
the orbitals: v — 1 (xz), v = 2 (yz), v = 3 (xy), v = 4 
(x 2 — y 2 ), and v = 5 (3z 2 — r 2 ). 



III. RESULTS 

In this section we present a systematic analysis of the 
magnetic beh aviour of the two-orbital model of Raghu et 
a/™ (section IIIA[), th e three-orbital model of Daghofer 
et alW^ (sect ion IIIB ), the four-orbital model of Yu et 
q/P^ (section III CI), a nd the five-orbital models of Kuroki 
et a/.™ (sectionW™ and Graser et a/. 48 (section JIBE). 



In all cases, we first determine the mean-field phase di- 
agram using the method outlined in section II A[ and 
then calculate the susceptibilities at several points in 
the PM state lying close to the boundary of the ordered 
phases. We emphasize that we do not seek to construct 
the true mean-field phase diagram of the orbital mod- 
els, but rather to determine the appropriateness of the 
mean-field ansatz O). 

All calculations (i.e. the phase diagram and the mean- 
field susceptibilities) were performed using a 400 x 400 k- 
point mesh. In order to better distinguish the total and 
the orbitally resolved susceptibilities from one another, 
we use a different color scheme in the density plots of 
these quantities. We also note that the total static sus- 
ceptibility is symmetric about the line q x = q y ; if peaks 
in x '"(q? oj = 0) are found off this line, for simplicity we 
will only mention the peak with q x > q y in the discussion. 



y ' \v ' vj 



0.5 

> ° 

-0.5 



1500 



i> O y C 



:v^ 



' J L l_ 



-1 -0.5 0.5 1 

(a) * (b) 






2 






^r~~^--^_____- 







_J\^_^^r^^___ 




___v y._ 


> 


-2 






\ / - 


i 


-4 






\/ 


w 


-6 






- 




-8 
1 n 






- 



(c) -r 



M 



FIG. 1: (a) The Fermi surface of the two-orbital model of 
Raghu et a/. 43 showing the dominant orbital contribution, (b) 
The orbital weight around the a and f3 Fermi surfaces as 
a function of the winding angle measured with respect to 
the y-axis, taken in the anti-clockwise direction. The orbital 
weight around the 7 Fermi surface is not shown as it has 
almost complete yz character, (c) The band structure along 
high-symmetry directions. 



A. The two-orbital model of Raghu et al. 

The two-orbital model of Raghu et al. 43 was one of 
the first attempts to model the pnictides from an orbital 
point of view, and al so the firs t of a number of simi- 
lar two-orbital modelsP 2 * 52 * 55 * 77 * As it only keeps the xz 
and yz orbit als, it may be regarded as a minimal orbital 
model of the system. Although it is likely too simple to 
give quantitative agreement with experiment, it is nev- 
ertheless of interest to us as it qualitatively captures the 
variation of the orbital character about the Fermi sur- 
face and hence provides important clues to the ordering 
mechanism. 

The Fermi surface and the b and s tructure of the two- 
orbital model is shown in figure [IJ 84 As we shall see, the 
Fermi surface is somewhat unlike those of the more so- 
phisticated models, specifically one of the xz/yz- derived 
hole pockets is located at (ir, tt) instead of at (0, 0). This 
implies rather different magnetic properties compared to 
these other models, as the nesting of the hole pocket at 
(7r, tt) with the electron pocket at (0, tt) allows both elec- 
tron pockets to participate on an almost equal footing in 
the (tt, 0) AFM state. 

The mean- field phase diagram of the two-orbital model 
is presented in figure ^J The (tt, 0) AFM state is stable 
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FIG. 2: The T-U mean- field phase diagram of the two-orbital 
model of Raghu et al.. 43 The location of the symbols indicate 
a co-ordinate at which the free energy was evaluated, and the 
corresponding minimizing state. The shaded region approx- 
imately indicates the extent of the (tt, 0) AFM phase. We 
show the static susceptibility at the points indicated by the 
star symbols in figure [3] 



above a critical interaction strength U c ~ 2.65 eV, with 
the critical temperature increasing linearly with U above 
this value. Our results are consistent with those of Ref.l50l 
for J = 0.1 £7. In a similar two-orbital model, 41 the mean- 
field ground state was found to be a superposition of 
Qx — (tt, 0) and Q y = (0, tt) stripe magnetic orders with 
mutually-perpendicular spin polarizations. Although our 
aim is to check the validity of the mean-field phase di- 
agram by examining the spin susceptibility just above 
the mean- field ordering temperature, we are not able to 
determine whether such a two-Q state or the expected 
(7r, 0) AFM phase is realized: in both cases we will find 
a peak in the susceptibility at (tt,0). 

In figure [3] we show the static susceptibility at the 
points indicated in figure [2J At both points the sus- 
ceptibility is sharply peaked at or very close to (tt, 0), 
reflecting the good nesting of the electron and hole pock- 
ets, and confirming the appropriateness of our mean-field 
ansatz. The absence of any other significant structure in 
the susceptibility implies that only scattering between 
the electron and hole pockets is important to the mag- 
netic ordering. Examining the susceptibility at ~ 150 K 
above the AFM transition temperature at a number of 
interaction strengths between U = 2.7 eV and 3eV, we 
find only minor quantitative differences in \ '"(q? °° — 0) 
within this range. At U < U c the magnetic response in- 
dicates weak incommensuration, with the susceptibility 
peaking at q ~ (0.957T, 0) [see figure[3la)]. An incommen- 
surate AFM phase is in fact realized at U — 2.6 eV and 
T = 100 K, as revealed by the observation of negative val- 
ues of x h ( ( l^ = 0) here. The unrestricted real-space 
Hartree-Fock of Ref. |41] may have missed the weakly- 
incommensurate phase because the considered cluster 
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FIG. 3: The static transverse spin susceptibility of the two-orbital model of Raghu et al~ai (a) U — 2.5 eV and T — 20 K and 
(b) U = 2.7 eV and T = 400 K. (c) The dominant contribution to the static susceptibility from the xz-xz term at U — 2.7 eV 
and T — 400 K. Note the different scales used in each panel. 



size was too small. 

The dominant contribution to the static transverse 
susceptibility is from the xz-xz and yz-yz terms. 
Xi^!(q^ = 0) is shown in figure J^c); X^ti^^ = °) can 
be obtained from this by interchanging q x and q y . We 
see that although there are definite peaks at both (tt, 0) 
and (0,7r), the latter peak is somewhat higher than the 
former. This implies that the yz [xz] orbital is slightly 
more important to (ir,Q) [(0,tt)] AFM order, consistent 
with the mean- field calculation in Ref. l50l 



B. The three-orbital model of Daghofer et al. 

The three-orbital model of Daghofer et al. 59 is the 
newest model studied in this work. In addition to the xz 
and yz orbitals, this model also includes the xy orbital. 
Assuming that the magnetic order is driven by nesting 
of electron and hole Fermi surfaces, it is reasonable to 
neglect the x 2 — y 2 and 3z 2 — r 2 orbitals as ab initio cal- 
culations predict that they have very little weight at the 
Fermi energyP^^ 

The electronic structure of the three-orbital model is 
summarized in figure |4j The model reproduces key fea- 
tures of the Fermi surface predicted by ab initio calcu- 
lations: there are two xz/yz-derived nearly circular hole 
pockets at the centre of the Brillouin zone, and a yz/xy- 
derived [xz/yz-derived] elliptical electron pocket centred 
at k = {tt, 0) [k = (0,7r)]. We note that the hole pock- 
ets have much greater effective mass than the electron 
pockets, in striking contrast to the two-orbital model dis- 
cussed above. 

The ground state phase diagram of the model as a func- 
tion of J and U was mapped in Ref. 59 using a mean-field 
ansatz that also allowed for stripe or staggered orbital or- 
der in addition to FM and AFM states. At moderate- to 
strong-coupling, orbital ordered states were obtained; in 
the weak-coupling regime and J = 0.25£7, however, the 
(7r, 0) AFM state without orbital order was found to be 
stable. The authors of Ref. [64] mapped the same phase 
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FIG. 4: (a) The Fermi surface of the three-orbital model of 
Daghofer et a/. 5 ^ showing the dominant orbital contribution, 
(b) The orbital weight around the labeled Fermi surfaces, (c) 
The band structure along high- symmetry directions. 



diagram using a more restricted mean- field ansatz, but 
their RPA calculations suggest that the (n, 0) AFM state 
is unstable towards incommensurate order. In figure [5j 
we show the phase diagram in the T-U plane. 

Our results for the spin susceptibility are in agreement 
with Ref. [64j and we therefore conclude that the phase 
diagram figure [5] does not reflect the actual mean-field 
behaviour of the model. At relatively weak coupling the 
dominant magnetic response occurs in a ring of almost 
uniform height centred at q = (0, 0) of radius ~ 0.5tt and 
0.4tt, see figure p[a) and (b) respectively. In the former 
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FIG. 5: The T-U mean- field phase diagram of the three- 
orbital model of Daghofer et aZ.P"P, The shaded region ap- 
proximately indicates the extent of the (tt, 0) AFM phase. 
We show the static susceptibility at the points indicated by 
the star symbols in figure [6] 



case, we can make out another, fainter ring of radius 
~ 0.657r; the radii of these rings correspond exactly to 
twice the average radii of the hole pockets, indicating 
that the magnetic response is dominated by scattering 
from one side of the hole Fermi surfaces to the other. The 
smaller radius of the ring in the U = 0.68 eV data reflects 
the reduction of the size of the heavy hole pockets by the 
increase of the chemical potential with temperature. In 
both of these figures, there is only a very small peak 
close to (7r,0); this peak is barely visible in the U = 
0.68 eV, T = 350 K data. The form of the susceptibility 
changes significantly at higher U: as shown in figure[6lc), 
two rather broad peaks develop at q ~ (0.757T, 0). The 
origin of this behaviour is unclear, as this wavevector 
does not seem to correlate with any Fermi surface or 
band structure feature. 

It is interesting to consider the consequence of shrink- 
ing the hole Fermi surfaces in this model, for example 
by electron-doping. We expect that this should also 
shrink the ring feature in the susceptibility, leading to 
significant weight close to q = 0, favouring FM order 
and triplet superconductivity. Such a situation might 
be relevant to the stochiometric superconductor LiFeAs, 
where ARPES reveals small heavy hole pockets which are 
poorly nested with electron Fermi surfaces.^ A model for 
this compound which captures these salient band struc- 
ture features was recently shown to display FM order and 
triplet superconductivity in the weak-coupling limit PS 
Intriguingly, unpublished NMR results on single crys- 
tals of LiFeAs report that there is no suppression of the 
Knight shift below T c for magnetic fields perpendicular 
to the crystal c axis, strongly suggesting triplet pairing.^ 

We note that there are several alternative three-orbital 
models which display an ^-dominated hole pocket at the 
M point P=H3J A similar analysis to that performed here 



has been carried out for such a three-orbital model, 62 
and the susceptibility was found to be peaked at an 
incommensurate wavevector close to (jr, 0). More re- 
cently, Zhou and Wang studied a three-orbital model of 
LaFeAsO using Gutzwiller mean-field theory restricted 
to two-site magnetic unit cells and found significant de- 
viations from the usual Hartree-Fock approach. 65 



The four-orbital model of Yu et al. 



The four-orbital model of Yu et a/. 54 includes all but 



y 2 or 



the 3z 2 — r 2 orbital. Although neither the a 

the 3z 2 — r 2 orbitals significantly contribute to the Fermi 

surface, the former is filled at lower doping levels. This 

suggests that the 3z 2 — r 2 orbital will be less relevant 

to the magnetism} 5 -^, justifying its exclusion from this 

model. 

As shown in figure [7| the structure of the Fermi sur- 
face is broadly similar to that of the three-orbital model, 
with elliptical electron Fermi pockets at the X points, 
and two nearly circular hole pockets with much greater 
effective mass at the T point. Note, however, that the 
electron pockets are rotated and have reversed dominant 
orbital content (i.e. xz «-»• yz), while each Fermi surface 
is smaller. These are likely relatively superficial differ- 
ences: the orientation of the electron pockets is correct 
in the (physical) two-Fe Brillouin zone, while the size of 
the Fermi pockets can be tuned through the tight-binding 
parameters. More seriously, however, the model contra- 
dicts some key predictions of ab initio calculations : the 
xy weight at the Fermi surface is heavily suppressed ) 71 * 72 * 
and only a single band contributes to the electron Fermi 
pocketsPEU 

Our mean- field phase diagram agrees with the T = K 
analysis of Yu et al: we find a low-temperature (it, 0) 
AFM state for 1.9eV < U < 2.65eV, while above a crit- 
ical U c ~ 2.65 eV the ground state is the (tt, tt) AFM 
which has a very high ordering temperature (> 3000 K). 
It is possible that the strong-coupling (tt, 0) AFM state 
found in Ref. |54] has lower free energy for U > 4 eV, but 
at these coupling strengths the AFM states are insulat- 
ing and hence irrelevant to the magnetic ordering in the 
pnictides. 

In figure [9] we show the static transverse susceptibil- 
ity at the points indicated in the phase diagram fig- 
ure[8J The \ow-U results [ngure^a) and (b)] are strongly 
reminiscent of the two-orbital model. For U = 1.8 eV 
and T = 20 K, just below the critical coupling strength 
of the (7r, 0) AFM state, \ '"(q?^ = 0) is peaked at 
q w (0.95tt,0). Closer to the (tt, 0) AFM phase there is 
evidence that an ordered incommensurate state is real- 
ized. At U = 2.1 eV and T = 400 K, we find a sharp 
peak exactly at q = (71", 0). These peaks are clearly de- 
rived from the very good nesting of the electron and hole 
Fermi surfaces. This conclusion is supported by examina- 
tion of the dominant contributions to the peak at (7r, 0) 
in figure p)Td-f): the xz-xz term has the largest value, 
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FIG. 6: The static transverse spin susceptibility of the three-orbital model of Daghofer et a/. 59 at (a) U ■ 
(b) U = 0.68 eV and T = 350 K, and (c) U = 0.88 eV and T = 980 K. 
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FIG. 7: (a) The Fermi surface of the four-orbital model of 
Yu et alr^* showing the dominant orbital contribution, (b) 
The orbital weight around the labeled Fermi surfaces in (a), 
(c) The band structure along high-symmetry directions. 



but the sum of the xz-yz and xz-xy terms is comparable. 
This is to be expected, as the xz orbital dominates the 
electron pocket at k = (tt,0) [see figure [7]; furthermore, 
the elliptical form of the electron pocket means that it is 
nested with the xz-dominated regions of the hole pock- 
ets, hence accounting for the dominance of Xit- There is 
nevertheless significant yz and xz orbital weight around 
the hole and electron pockets, respectively, accounting for 
the importance of the inter-orbital susceptibilities \it 

and xi~,3~- 

Consistent with the mean-field phase diagram, a broad 



peak develops in x '"(q?^ = 0) at q = (71", n) with in- 
creasing U; close to the (7r,7r) AFM state it dominates 
the magnetic response, see figure [9jc) . Although there 
is low-energy scattering between the two electron pock- 
ets with this wave vector, the large magnitudes of the 
mean-field staggered magnetic potentials r\ v [([6|] sug- 
gests that the (7r,7r) AFM state cannot be understood 
only in terms of Fermi surface physics. For example, at 
U = 2.7 eV and 3eV we have 0.53 eV> \r] u \ > 0.43 eV 
and 1.23 eV> \r] u \ > 0.93 eV, respectively. The former is 
comparable to the minimum energy (at the X points) of 
the band lying entirely above the Fermi surface; the lat- 
ter is close to the difference between the hole bands and 
the Fermi energy at the M point. Since the r\ v are the 
characteristic energies for the reconstruction of the elec- 
tronic structure in the AFM state, it is therefore likely 
that these regions participate in the (tt, tt) AFM state, 
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FIG. 9: The static transverse spin susceptibility of the four-orbital model of Yu et atf^at (a) U — 1.8 eV and T = 20 K, (b) 
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accounting for its much lower free energy compared to 
the (tt, 0) AFM phaseP 



D. The five-orbital model of Kuroki et al. 

The five-orbital model of Kuroki et alW^ includes all 
five Fe d-orbitals. It was constructed by fitting to ah ini- 
tio calculations for LaFeAsO, and it was one of the first 
orbital-based models proposed for the iron pnictides. The 
behaviour of this model, and several of its close relatives, 
has been extensively studied j^^ 49 ' 51 ' 56 ' 57 ' 60 ' the mag- 
netic properties are particularly well understood, with 
RPA calculations confirming that (7r, 0) AFM order is 
realized in the undoped model P*^ The extent of this 
state in the phase diagram at zero doping, and the possi- 
ble competition with other magnetic phases, nevertheless 
remains largely unexplored. 

This model was originally constructed in an orbital ba- 
sis where the x and y axes are rotated by 45° to the Fe 
lattice. In the other models considered here, however, the 
Fe lattice defines the coordinate axes for the orbital basis. 
To allow direct comparison with these other models, we 
therefore show the Fermi surface of the five-orbital model 
of Kuroki et al. in the usual basis in figure [lO^a). Al- 
though the Fermi surface has the same basic form and or- 
bital composition as the three-orbital model, the electron 
pockets are almost circular, and there is an additional 



^-dominated hole pocket at the M point. Furthermore, 
from the band structure figure [lQ^c) we see that the hole 
pockets at the T point have similar Fermi velocities to 
the electron pockets, which themselves have much higher 
Fermi velocity along the T— X line than in the X— M di- 
rection. The latter accounts for the anisotropic spin- wave 
dispersion found in Ref. [571 



In constructing our mean- field phase diagram figure [TT 
we extended the ansatz pi) to also compare the free en- 
ergy of the ordered states within the rotated orbital basis 
used by Kuroki et al. to those in the usual orbital ba- 
sis. We find that the mean- field solution in the usual 
orbital basis always has the same [PM and (jr, ir) AFM] 
or slightly lower [(tt, 0) AFM] free energy than the mean- 
field solution in the rotated basis. Although our results 
are not directly comparable with those of Kaneshita et 
a/.^who used the rotated orbital basis and allowed for 
inter-orbital mean fields, our phase diagram is neverthe- 
less in good qualitative agreement ! 4T | 49 1 for U c > 0.85 eV 
we find a (7r,0) AFM state whose critical temperature 
rapidly grows with increasing U. At U > 1.15 eV and 
T > 1000 K, however, the (7r,7r) AFM state is stable. 

In figure [12| we present the results for the static trans- 
verse susceptibility at the points marked in figure 11 
For U < 1 eV, the results are rather similar to those al- 
ready seen in the two- and four-orbital models. In the 
low-temperature PM state [figure [l2^ a)], the susceptibil- 
ity takes a maximum in an arc of radius ~ 0.057T about 
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FIG. 10: (a) The Fermi surface of the five-orbital model of 
Kuroki et a/= showing the dominant orbital contribution, 
(b) The orbital weight around the a, f3 and 7 Fermi surfaces 
in (a); the orbital weight around the 5 Fermi surface is not 
shown as it has almost complete xy character, (c) The band 
structure along high-symmetry directions. 
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FIG. 11: The T-U mean- field phase diagram of the five- 
orbital model of Kuroki et a/.. 40 The lightly-shaded region 
approximately indicates the extent of the (71", 0) AFM phase, 
while the darkly-shaded region gives the extent of the (7r, 7r) 
AFM phase. We show the static susceptibility at the points 
indicated by the star symbols in figure [12] 



(7r, 0); at slightly higher U = 0.85 eV we find evidence 
of an incommensurate AFM state with ordering vector 
(0.9257T, 0.0257r). For the points lying directly above the 
mean-field (71", 0) AFM state, however, the susceptibil- 
ity takes a maximum exactly at q = (71", 0). This peak 
is nevertheless much broader than in the two- or four- 
orbital models, indicating a larger distribution of nesting 
vectors. 

The dominant contributions to the peak at q = (71", 0) 
in the U = 0.9 eV, T = 300 K static susceptibility [fig- 
12 'b)] are shown in figure fl3] The most important 



ure 



terms are the yz-yz [figure [l3[a)j, yz-xy, [figure [l3|b)] 



and xy-xy [figure [l3[c)] susceptibilities, which together 
constitute « 50% of the peak height in the total sus- 
ceptibility; for comparison, in the four-orbital model the 
dominant susceptibilities shown in figure (9J^d-f) are re- 
sponsible for « 90% of the (71", 0) peak. The dominant 
terms in Kuroki et a/.'s five-orbital model are consistent 
with nesting of the outer hole pocket with the electron 
pockets being the main driver of the magnetism, but with 
the nesting of the hole pocket at the M point with the 
electron pockets also playing a substantial role. This 
is supported by the reconstructed Fermi surface shown 
in Ref. HTlfor the low-moment case. 

As U is increased to 1.05 eV [figure (12^ d)] the nesting 
picture of the magnetism begins to break down, with the 
appearance of a definite (but still subdominant) peak at 
q = (7r,0.57r). This feature completely dominates the 
response at U = 1.1 eV, indicating that the mean-field 
phase diagram is unreliable for higher coupling strengths. 
Peaks at these wavevectors have previously been ob- 
served pUaSI where they were explained as due to scatter- 
ing between the electron Fermi pockets. This appears to 
be inconsistent with the observation that these peaks are 
most pronounced at high temperatures > 1000 K, which 
instead suggests that band structure features away from 
the Fermi surface are responsible. In particular, we note 
that there is a region of flat bands at —0.2 eV near the M 
point, which is predo minantl y derived from the almost- 
filled 3z 2 - r 2 orbital-SSEHSD These states are shifted to 
higher energies by the Hartree terms, and hence might 
play a major role in the magnetic response at high tem- 



peratures. This will be discussed in detail in section IV B 



E. The five orbital model of Graser et al. 

The five-orbital model of Graser et aZ.™ is a frequently 
studied alternative to Kuroki et a/.'s model. Graser et 
a/.'s model was also constructed by fitting to ah initio re- 
sults; as can be seen from the electronic structure shown 
in figure [l4j both five-orbital models share many similar- 
ities. The key differences in the Fermi surface are that 
the electron pockets in this model are clearly elliptical, 
while the hole pocket at the M point is almost vanish- 
ingly small. 87 The former is in much better agreement 
with current ARPES and ah initio results; the latter 
refinement is also in reasonable agreement with ah ini- 
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FIG. 12: The total static transverse spin susceptibility of the five-orbital model of Kuroki et a/. 40 at (a) U — 0.8 eV and 
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FIG. 13: The dominant contributions to the (71", 0)-peak at U = 0.9 eV and T — 300 K in the five-orbital model of Kuroki et 
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tio calculations, but the numb er of hole p ockets remains 
somewhat unclear in ARPESP^ 71 ' 72 ' 88 ' The flat band 
region near the M point is also closer to the Fermi surface 
than in Kuroki et a/.'s model; as such, we might expect 
it to play a larger role in the magnetic response. 

The magnetic properties of Graser et a/.'s five-orbital 
model has only recently attracted attention. In Ref. [61] 
it was shown that the (tt, 0) AFM state was realized at 
T = K for U > 1.23 eV, but the authors did not consider 
competition with other magnetic states. The existence of 
a (71", 0) AFM ground state with realistic ordered moment 
and Fermi surface was confirmed in Ref. 64 for a rather 



narrow range of U. For much larger values of U it has 
been shown that Hartree-Fock and Gutzwiller theories 
give quite different results. 68 In figure 15 we show the 



phase diagram in the T-U plane: a (tt, 0) AFM state with 
low critical temperature < 500 K is realized for 1.25 eV< 
U <1.35eV, but the phase diagram is clearly dominated 
by a (7r,7r) AFM state with high critical temperature 
> 2000 K. For most of the U range where the (tt, 0) AFM 
state is stable at T = OK, we find that the (tt, tt) AFM 
state is in fact realized at higher T, with the PM phase 
separating them at lower values of U. 



The results for the static susceptibility [figure 16 in 
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AFM phase. We show the static susceptibility at the points 
indicated by the star symbols in figure [16] 



dicate that the mean-field phase diagram is only par- 
tially correct. The evolution of x ^(o.^ — 0) U P to 
the edge of the (71", 0) AFM phase is similar to that 
in Kuroki et a/.'s model: as U is increased, a definite 
tendency to weakly-incommensurate order [with a peak 
at q « (0.957r,0.0257r)] gives way to broad peaks at 
q = (71", 0). We have confirmed an instability towards in- 
commensurate order below at least 100 K at U — 1.2 eV. 
We also note a ring-like feature centred at q = (ir, tt) 
in figure [l6^a); by U = 1.23 eV this has apparently 
evolved into a broad peak at q = (7r,0.67r). Although 
this is similar to the additional peaks seen in the high- 
temperature results for Kuroki et a/.'s model, they ap- 
pear here at much lower temperatures. 

The peak at (71- , 0.67r) displays a peculiar temperature- 
dependence. To see this, in figure 16 ^d) and (e) we show 
the static susceptibility at U — 1.25 eV for T = 500 K 
and 900 K, respectively. In the former case, the peaks at 
both (7T, 0) and (n, 0.6ir) are distinctly visible. We have 
verified that as the temperature is lowered to ~ 450 K the 
(7T, 0) peak diverges while the (tt, 0.6tt) peak is somewhat 
suppressed, and so the (tt, 0) AFM state is realized at 
low temperatures in agreement with the mean-field phase 
diagram. As the temperature is raised to T = 900 K, 
however, the peak at (tt, 0) disappears while the (tt, 0.6tt) 
peak almost doubles in height, revealing that the system 
is close to an instability towards a (7r, 0.67r) AFM state. 
This also suggests that the mechanism responsible for the 
(tt, 0.67r) peak is not due to scattering between different 
Fermi pockets, but rather due to the electronic structure 
away from the Fermi surface. In figure 16 'f ) we show the 
static susceptibility at slightly higher U and T compared 
to figure [l6^d): the (7r,0.67r) peak completely dominates 
the response, and the system is again close to ordering 
at this wavevector. Along with fig ure [l6^ e), this suggests 
that the (71", tt) AFM state in figure [l5fsnould be replaced 
by a (tt, 0.6tt) AFM state; the region of stable (tt, 0) AFM 
order will also likely shrink. 

To investigate the origin of the (jr, 0.6tt) peak we have 
examined the dominant contributions to \ V - Focusing 
upon the U = 1.25 eV, T = 500 K case, we find that the 
susceptibilities involving the xy-oibital are most impor- 
tant, contributing « 50% of the peak height; the largest 
contributions come from the xy-xy term [figure (17^ a )], 
the yz-xy term [figure |l7|fo)], the xz-xy term [figur e |17^ b) 
with q x o q y ], and the (x 2 — y 2 )-xy term [figure |l7[c)]. 
The dominant role of the xy orbital is consistent with 
the results of Ref. [51] and Ref. [56 for similar orbital 
models, but it is difficult to reconcile the temperature- 
dependence observed here with the author's interpreta- 
tion that the (71", 0.67r) peak originates from scattering 
between the electron pockets. An alternative possibility 
is scattering between the hole pockets at the Y point and 
the ^-derived hole pocket at the M point, although the 
required wavevector is slightly too large. We will discuss 
this matter further in section HVBI 
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FIG. 16: The total static transverse spin susceptibility of the five-orbital model of Graser et al. 48 at (a) U = 1.15 eV and 
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IV. DISCUSSION 



To summarize our main results, we found that the two- 
orbital, 43 four-orbital, 54 and Kuroki et a/.'s five-orbital 40 
models display a rather robust (ir, 0) AFM phase in their 
zero-doping phase diagrams. In both the four- and five- 
orbital models, this AFM state competes with a strong- 
coupling AFM state with different ordering vector. In 
contrast, the three-orbital 5 -^ and Graser et a/.'s five- 
orbital 48 models show at most a weak tendency towards 
(7r, 0) AFM order: in the former it is likely not present 
anywhere in the considered phase diagram, while in the 



latter it is out-competed over much of the phase diagram 
by a much more stable (7r, 0.67r) AFM state. 

Our results support the scenario of the (jr, 0) AFM 
state originating from nesting of the electron and hole 
Fermi pockets. This is a crucial ingredient for the mag- 



netism: as discussed below in section IV A, the stabil- 



ity of the (7r, 0) AFM state is positively correlated with 
the degree of the nesting. The similarity of the elec- 
tronic structure of the two five-orbital models, yet the 
very large differences in their magnetic phase diagrams, 
is also of special note. In particular, this sheds light on 
the origin of the (7r, 0.57r)-(7r,0.67r) AFM state realized 
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in these models at strong coupling or high temperatures. 



In section |IVB| we argue that subtle differences in the 
states derived from the 3z 2 — r 2 orbital play a significant 
role in stabilizing this state. Finally, in section IV C we 



discuss the implications of our results for models of the 
magnetic order in the iron pnictides. 



A. Nesting 

The (7r, 0) AFM order was found to be most robust 
in models where the quality of the nesting between the 
electron and hole pockets is high. For perfect nesting the 
dispersion ek satisfies 



ek 



~ e k+Q 



(17) 



for some nesting vector Q. In such a scenario, not only 
do the Fermi surfaces exactly overlap upon performing 
the shift by Q, but also the holes and electrons have the 
same effective mass. In realistic models such as those 



considered here the condition (17) is only approximately 



fulfilled. For a nesting picture to still be relevant, we re- 
quire that some segments of the electron and hole Fermi 
surfaces should overlap, and that the effective mass dif- 
ference between them should not be too great. We note 
that it is sometimes the case that while the condition (17) 



is satisfied approximately for the entire Fermi surface for 
the vector Q, a better match between some Fermi sur- 
face segments (but worse for others) may be achieved 
for a slightly different vector Q « Q. At weak coupling 
strengths, where the effective staggered magnetic poten- 
tial is small, more of the Fermi surface may be gapped 
at ordering vector Q than at Q. 

In figure [IHJ we show the Fermi surface of each of the 
models superimposed with the Fermi surface translated 
by Q = (tTjO). It is immediately clear that the three- 
orbital model has the worst nesting properties: not only 
is there very poor overlap between the hole Fermi surface 
and the translated electron Fermi surface, but also the 
electron pocket has much lower effective mass than either 
of the hole pockets. It is therefore not surprising that this 
is the only model that fails to show the required (7r, 0) 
AFM order. 

Although both electron pockets participate in the 
(7r,0) AFM state in the two-orbital model, the nesting 
here is nevertheless rather similar to that in the four- 
orbital model. In these two models, a segment near the 
minor axis of the shifted elliptical electron pocket over- 
laps with the small hole pocket, while a segment near the 
major axis overlaps with the large hole pocket [see fig- 
ure fl8[a) and (c)]. There is an excellent match between 
the dominant orbital character of the original and the 
shifted Fermi surfaces, which should enhance the AFM 
fluctuations as the intra-orbital Coulomb repulsion is 
much larger than the inter-orbital Hund's rule coupling. 

The nesting of the electron pockets with the outer hole 
pocket at the T point in Kuroki et aZ.'s five-orbital model 



is the best considered here, with almost complete overlap 
with the shifted Fermi surfaces. The orbital-resolved sus- 
ceptibilities suggests that this is most important to the 
magnetic instability, although the poorer nesting of the 
electron pockets with the hole pocket at the M point is 
likely responsible for broadening the peak at (7r,0). Un- 
like the two- and four-orbital models, therefore, there is 
a relatively poor match of the dominant orbital character 
of the nested Fermi surfaces. This is reflected in the sus- 
ceptibilities: ^2 U Xut contributes w 30% of the height of 
the (7r, 0) peak in Kuroki et a/.'s model, while in the four- 
orbital model the susceptibility X;Ti~ * s a l° ne responsible 
for « 50% of the peak height. In contrast, despite the 
near-identical orbital content of the Fermi surfaces, the 
instability towards (jr, 0) AFM order in the five-orbital 
model of Graser et al. is much weaker. This can be 
explained by the absence of direct Fermi-surface overlap 
between the electron and the hole pockets at the T point, 
and also the much weaker nesting between electron pock- 
ets and the small hole pocket at the M point. The stark 
differences between the two five-orbital models highlights 
the sensitivity of the (ir, 0) AFM state to subtle details 
of the electronic structure. 

Our results suggest that a qualitatively-correct under- 
standing of the AFM order in the pnictides can be devel- 
oped without reference to the orbital structure, i.e. they 
support the excitonic theories of the magnetism. This 
is not unexpected, as the (-zr, 0) AFM order is realized 
at relatively weak coupling strengths, where the exci- 
tonic model arises naturally as a low-energy effective the- 
orv p9|33 | ^he ra ther strong competition with other AFM 
states, however, is not anticipated by such excitonic mod- 
els. 



B. Hartree shifts 

The inclusion of Hartree shifts in our mean-field de- 
coupling scheme requires justification: models where the 
tight-binding band structure is constructed by fitting to 
a b initio results should already account for the Hartree 
shifts, so it is then only necessary to include the stag- 
gered magnetic potentials in Q. It is nevertheless likely 
that such methods underestimate the effect of the correla- 
tions; indeed, comparison with ARPES reveals that the 
bands obtained within the local density approximation 
are subject to significant energy shifts. 12 Including the 
Hartree shifts is the simplest way to test the robustness 
of the AFM order towards such modifications of the band 
structure. In Ref. [56] it was reported that the Hartree 
shifts can lead to dramatic changes in the Fermi surface 
of Kuroki et a/.'s five-orbital model. Ikeda and co-workers 
similarly found unphysical renormalization of the Fermi 
surface by the uo = self-energy corrections calculated 
within the fluctuation-exchange approximation. ^ * 51 * 60 * 

In order to evaluate the effect of the Hartree shifts, 
we have calculated the renormalized band structures for 
each model in a fictitious T = K PM state at a coupling 
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strength where an AFM mean-field state is realized. This 
is a good estimate of the band structure in the stable T > 
K PM state as there is only rather small temperature- 
dependence of the orbital occupancies. The degeneracy 
of the xz and yz orbitals in the two-orbital model means 
that there is no change in the band structure; for the 
three- and four-orbital models we find only very slight 
changes in the band structure far from the Fermi surface 
(> 0.5 eV and leV, respectively). As such, the inclusion 
of the Hartree shifts should be irrelevant to the weak- 
coupling magnetic response of these models. 

For the five-orbital models we again find very minor 
changes in the Fermi surface, limited to the shrinking 
(Kuroki et al.) or removal (Graser et al.) of the xy- 
derived hole pocket at the M point, which is due to 
the positive Hartree shift of the other orbitals relative 
to the (least-occupied) xy orbital. As shown in fig- 
ure Il9j a) and (c), however, the flat-band region near the 
M point and below the Fermi energy undergoes an al- 
most uniform shift to higher energies. These states are 
mostly derived from the 3z 2 — r 2 orbital, which has the 
highest occupation and hence has the largest Hartree 
shift li§ll 8 l 51 l 60 l The shift to higher energies is particu- 



larly dramatic in Graser et a/.'s model, where the fiat 
band region is shifted from « — 0.15eV to « — 0.07eV; 
in Kuroki et a/.'s model the equivalent feature is shifted 
from « -0.21 eV to « -0.17 eV. We note that the 
Hartree shifts of the orbitals are in qualitative agreement 
with the more advanced analysis of Ikeda et al. j 46 | 51 | 6Q | al- 
though the change in the band structure is less severe in 
our results. 

The shift of the (3z 2 — r 2 )-derived states closer to the 
Fermi surface suggests that they might play a role in 
the high-T and U magnetic behaviour of the system, 
where the {tt, 0) AFM state competes with an enigmatic 
(7r,0.57r)-(7r,0.67r) AF M state . This is contradicted by 
the result from section HIE that the peak at (7r,0.67r) 
is mostly due to the xy orbital, with the 3z 2 — r 2 or- 
bital making a relatively modest contribution. To isolate 
the effect of the 3z 2 — r 2 orbital, we have therefore re- 
calculated the static susceptibility without including the 



Hartree shifts, x + (q,cc; = 0). In figure 19 'b) and (d) 



we show the ratio of the static susceptibility with the 
Hartree shift to that without for the cases figure [l2^d) 
(U = 1.05 eV, T = 1000 K, Kuroki et a/.'s model) and fig- 
ure [l6jd) (U = 1.25 eV, T = 500K, Graser et a/.'s 
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model). In both cases, we see that the peak at (7r,0.57r)- 
(7r, 0.67r) is significantly enhanced by the inclusion of the 
Hartree shifts, especially so for Graser et a/.'s model. In- 
cluding the Hartree shifts also leads to a reduction of the 
height of the peak at (jr, 0), probably due to the shrinking 
of the xy pocket, highlighting the sensitive dependence 
of the (7r, 0) AFM state upon fine details of the nesting. 

As the rest of the band structure undergoes relatively 
little change upon the inclusion of the Hartree shifts, the 
results figure 19 'b) and (d) indicate that the shift of the 



3z — r -derived flat band near the M point is heavily 
involved in stabilizing the (jr, 0.57r)-(7r, 0.6tt) AFM state. 
This is paradoxical: although the 3z 2 — r 2 orbitals are 
clearly important to the appearance of the (7r,0.57r)- 
(7r, 0.67r) AFM state, other orbitals dominate the mag- 
netic response at the ordering vector. Furthermore, the 
inclusion of the Hartree shifts enhances all the orbitally 
resolved susceptibilities %~+ by a similar factor. For ex- 
ample, in the U = 1.25 eV, T = 500 K results for Graser 
et a/.'s model, the ratio 2.3 > XuX /Xat ^ 1-8 for all ji 



and ^, although it is maximal when the 3z 2 — r 2 orbital 
is involved. This appears to indicate that the flat band 
near the M point plays an indirect role in the magnetic 
ordering. We speculate that this involves the strong peak 
in the density of states associated with this band, which 
feeds into the response of the other bands due to the 
orbital mixing. 



The sensitivity of the (n, 0.57r)-(7r, 0.67r) peak to the 
Hartree shift implies a strong dependency upon the 
choice of the ratio J/U. For J/U = 0.25 considered 
here, we have positive effective Hartree shifts e v = 
Q.2bU(n v — min{n M }). Had we chosen J/U = 0.2, how- 
ever, the Hartree shifts would vanish; for J/U < 0.2 we 
find negative Hartree shifts, and the 3z 2 —r 2 orbital would 
therefore be shifted away from the Fermi energy, likely 
suppressing the (7r, 0.57r)-(7r, 0.6tt) AFM order. 
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C. Implications for models of the iron pnictides 

One of the aims of this work was to identify an ideal 
low-energy orbital model of the pnictides, i.e. a model 
possessing Fermi surfaces in good quantitative agreement 
with ARPES and ab initio results, and also displaying 
robust (7r, 0) order. Unfortunately, none of the models 
studied here satisfy both criteria. The models which 
display robust (71", 0) AFM order have rather unrealis- 
tic elements in their electronic structure: the two-orbital 
model has a xz/yz- derived hole pocket at the M point 
and also neglects the large xy weight at the Fermi en- 
ergy, in the four-orbital model the xy weight is also 
unrealistically small at the Fermi surface and the elec- 
tron pockets do not involve a band crossing, while the 
electron pockets in Kuroki et aVs five-orbital model are 
not elliptical. In contrast, the three-orbital model qual- 
itatively captures the predicted Fermi surface topology 
and dominant orbital contributions, while Graser et a/.'s 
five-orbital model is in good quantitative agreement, but 
these two models show the weakest tendencies towards 
(tt, 0) AFM order. 

The sensitivity of the (71- , 0) AFM state to small details 
of the band structure raises an important question about 
the magnetic state in the pnictides: if electron-hole nest- 
ing is indeed responsible for this state, why is it realized 
in such a large range of compounds, each with subtly 
different Fermi surfaces'?^ A possible answer is that the 
coupling to the lattice degrees of freedom allows the sys- 
tem to fine-tune its Fermi surface. Indeed, the magnetic 
properties of these compounds have been shown to be 
strongly coupled to details of the lattice structure: the 
lattice constants and pnictogen height above the Fe plane 
obtained in ab initio calculations are strongly affected 
by the iron moment) 89 * 90 * and varying the Fe-A s-Fe bond 
angles can significantly alter the Fermi surfaced 6 * 77 * 78 * 91 * 
It is certainly suggestive that Kuroki et a/.'s model was 
fitted to ab initio results using experimental crystal pa- 
rameters, whereas Graser et a/.'s is based on structurally- 
relaxed calculations. 92 

Moreover, a structural phase transition from a high- 
T tetragonal to low-T orthorhombic phase occurs either 
just above (i?FeAsO compounds) or at (AFe2As2 com- 
pounds) the AFM transition temperature.™^ A scenario 
for a joint magnetic-structural transition in the pnictides 
was proposed within the context of a simple excitonic 
model in Ref. [93| the lattice distortion improves the 
nesting condition sufficiently to allow the appearance of 
an AFM state, which in turn stabilizes the distortion. 
This was verified for a more realistic four-band excitonic 
model in Ref. 32j Further work to clarify the role of the 
orthorhombic distortion is clearly necessary. It would be 
particularly interesting to examine this question within 
an orbital model, as the directional wave functions of the 
orbitals naturally couple to the lattice. The possibility 
that spontaneous orbital ordering drives the structural 
transition is of particular interest)- 94 * 95 * Of course, another 
possible explanation for the stability of the (jr, 0) AFM 



order involves strong-coupling physics beyond the nest- 
ing scenario! 65 * 68 * We note that recent calculations using 
dynamical mean-field theory suggest that the correlation 
strength in these materials has been underestimated. 96 
Further research is clearly needed to understand the com- 
plex magnetic physics of the pnictides. 



V. SUMMARY 

In this work we have presented a systematic weak- 
coupling investigation into the appearance of the ex- 
pected (tt, 0) AFM order in five different orbital models 
of the iron pnictide parent compounds. The mean-field 
phase diagrams of each model were first determined, and 
then the static spin susceptibility was calculated within 
RPA in the PM state close to the boundaries of the mag- 
netic phases. The highest peaks in the spin susceptibility 
reveal the actual ordering vector of the nearby mean-field 
state. The origin of these peaks was studied by examin- 
ing the dominant orbitally resolved contributions to the 
total spin susceptibility. This procedure allowed an un- 
biased assessment of the magnetic-ordering properties of 
each model. 

The studied models display an unexpectedly rich range 
of magnetic behaviour. Only four of the models were 
found to have an instability to an AFM state with 
the required ordering vector, although the robustness 
of this phase varies greatly between them. In each of 
these models we also uncovered evidence of a weakly- 
incommensurate AFM phase with ordering vector close 
to (7r,0) at coupling strengths below the minimum re- 
quired for a stable (71", 0) state. The (n, 0) AFM state 
originates from the nesting of the electron and hole Fermi 
pockets. The quality of the nesting appears to be of pri- 
mary relevance for the stability of the (n, 0) AFM order, 
while the matching of the orbital character of the nested 
Fermi surfaces is of lesser importance. In the most real- 
istic five-orbital models the (tt, 0) AFM state was found 
to compete with a strong-coupling incommensurate mag- 
netic state. We have argued that the origin of this phase 
involves electronic structure below the Fermi energy. The 
implications of our results for orbital models of the pnic- 
tides was discussed, and we propose that the apparent 
sensitivity of the (7r, 0) AFM order to fine details of the 
low-energy electronic structure indicates that lattice de- 
grees of freedom should be included in a theoretical de- 
scription of the magnetism of the pnictides. 
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